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1 Introduction 

In this paper we provide an analytical study of the theory of multi-valued and 
fuzzy cellular automata where the fuzziness appears as the result of the appli- 
cation of an underlying multi-valued or continuous logic as opposed to standard 
logic as used conventionally, e.g., |25j . Thus, we consider boolean elementary cel- 
lular automata (EGA) and in particular those CA considered by Wolfram |^ 
as a starting point. For an excellent review of classical CA theory and current 
problems within see [11]. Using the disjunctive normal form of any one of the 255 
ECA's so defined, or rules as we call them sometimes, we modify the underlying 
logic structure and redefine the ECA within the framework of this new logic. 
This then defines a CA under a new logic: For example, the classical ECA [25] . 
are CA under boolean (two-valued) logic, the fuzzy cellular automata (FCA) de- 
fined in [3] and considered in [3], [7], [S], [T5], [IS], [H], [H], are continuous CA 
under a special fuzzy logic (see Sections 3-4 herein). This latter logic has been 
used recently to generalize the game of Life, see e.g., PT , and is used to inves- 
tigate the result of perturbations, noisy sources, computation errors, mutations, 
on the evolution of boolean ECA, [8J. In addition, since all electromagnetic sig- 
nals travel at finite speeds, there is a point at which a boolean switch is neither 
ON nor OFF. Although this time interval is very short, this can be considered 
a fuzzy event, that is, an event where the switch is neither on nor off and is 
therefore in a countable number of intermediate states. As miniaturization and 
nanoscale technology becomes the norm, we suspect that such fuzzy or continu- 
ous states will play an important role in the development of future technologies, 
where quantum mechanical laws will intermingle with fuzzy states continuously 
and may even coincide. 

The idea here is to study the role of the logical structure on the dynam- 
ics of the new automaton and to determine whether or not complexity exists; 
in particular, the main question asked is: " Are the dynamics of the resulting 
space-time diagrams chaotic?" This is, of course, a difficult question to answer in 
this generality so we will target an important CA namely, ECA 110, and study 
its space-time diagram under a change of logic from boolean to probabilistic. 
Although it is somewhat surprising that FCA 110 is incapable of chaotic dy- 
namics (see [15]) as opposed to its boolean counterpart [25], it is definitely odd 
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that the present study of a special probabihstic CA known as PCA 110 also ap- 
pears to lean in this direction. One would expect that of all the continuous logics 
available, an underlying probabilistic logic (see Example 2) should yield chaotic 
space-time diagrams for the corresponding cellular automaton to boolean rule 
110, however, such a premise is not supported by theoretical arguments. The 
ultimate point of all this is to determine whether such non-chaotic evolution can 
shed some light into the sometimes chaotic evolution of the same automaton, 
under boolean logic, as the seed value approaches the number 1 from the left. 
We suspect that a necessary condition for the existence of chaotic space-time 
diagrams for a general fuzzy rule is that the underlying logic be discontinuous 
as a function of the two variables, in the sense of mathematical analysis. 

We will see that if the seed a, in a homogeneous background of zeros, is 
chosen sufficiently close to but less than 1, then the space-time diagrams (that 
is the space vs. time evolution of the automaton) of both boolean and fuzzy 
rule 110 (denoted by FCA 110) up to time T will be identical. Thus, one may 
expect that the asymptotic evolution of the fuzzy space-time diagram will have 
a bearing on the distribution of O's and I's in the corresponding asymptotic 
boolean evolution. More precisely, we can say that given any row number, or 
time, T > 0, there exists a seed value a < 1 such that the evolution of EGA 
110 under a specific fuzzy logic (dubbed CFMS logic, see below) up to row T 
is identical to its boolean evolution up to row T (provided we color all cells 
having values greater than one-half as black and all others as white) while for all 
subsequent rows the evolution is non-chaotic whereby various limits exist and 
are given by quotients of linear combinations of Fibonacci numbers, see |15j for 
specific details and the sequel for examples. We know that under CFMS logic all 
but 6 of the 255 EGA admit non-chaotic behavior for any finite random initial 
string [T7j, [TS]. A similar statement can be formulated for probabilistic EGA 110 
(denoted by PGA 110) but the precise result is conjectural at this time due to 
the massive amount of technical details that arise in its investigation (see the last 
section). We suspect that those GA that lead to truly random behavior (whether 
classical EGA or fuzzy) are those whose space-time evolution from one seed gives 
limits that are one-half, independently of the seed value. For example, EGA 18 
and its fuzzy counterpart, FGA 18 (see below) would enter into this realm. We 
can only surmise, and cannot prove completely, that PGA 110 is non-chaotic 
as well on finite initial strings. We leave it as an open question the problem of 
determining whether there is non-chaoticity (lack of complexity) in EGA under 
the other multi-valued logics mentioned here (Lukasiewicz, etc.) 

2 Logical structures and cellular automata 

Basically, a logic is multi-valued if it takes more than two values (e.g., more than 
just True or False, 1 or 0, etc.). Thus, by definition, multi- valued logics allow for 
truth- values that are more than just true or false. In a sense they allow values in 
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between, which is where the fuzziness concept comes in. For example, in classical 
propositional logic the law of the excluded middle states that for any proposition 
P, the statement {P V ^P) (P or not P) is always true. Such a statement allows 
for "proofs by contradiction" which is a technique that cannot be used generally 
for multi-valued logics. Since its origins in the 1920's there has been considerable 
work in multi-valued logics, in the past 20 years much of it resulting from the 
application of fuzzy techniques to real- world problems (cf., e.g., |22], |5]). For 
generalities concerning fuzzy automata we refer the reader to [[B], p. 199] while we 
cite |23j for details concerning fuzzy logics. Our aim here though is to entertain 
the use of general many- valued or continuous logics as well as general CA and to 
single out one or two of them in order to show that some fuzzy logics can give rise 
to non-chaotic behavior using the corresponding EGA. For various applications 
of many- valued logics to hardware design and artificial inteUigence, see [TD], [3], 

Whenever a logic takes on exactly three values it is called trivalent [|13j. 
p. 192]. The three values may take the form True, False, and something in between 
or we can also use numbers such as 1, 0, and 1/2 (cf., [H]). A logic that takes 
on finitely many values is called an n-valued logic while an infinitely- valued logic 
is one that takes on infinitely many values. Specifically, a continuous logic is 
one that takes on a continuum of values. In particular, in this article we will 
always assume that our logics have their values in the closed interval [0, 1] of 
real numbers. 

Denoting the two-point set (or boolean space) consisting of the numbers 
and 1 by {0, 1}, classical propositional (or standard) logic gives the truth tables 
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the truth table of ^P being found easily because of the law of the excluded 
middle. Recall that P Q is defined classically as (^-P) V Q. Prior to un- 
dertaking any constructions we illustrate by means of examples some of the 
multi- valued/continuous logics we are considering in this work. 

Example 1 Some common multi-valued logics and their formulations can be 
seen in the following two tables: The Lukasiewicz p-logic is defined for all p > 1, 
the case p = 1 being the most commonly used one. Since these logics have their 
values in the closed interval [0, 1], x,y G [0,1] therein. 

Remark 1. The logic described in the table below as CFMS logic seems to appear 
first in the paper [5]. Observe that the functions defined by on [0, 1] and xf\y, 
X V y on [0, 1]^ map back into [0, 1] for each logic presented. However, each of 
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Lukasiewicz p-logic 


Zadeh logic 


CFMS logic 




1-x 


l~x 


1 - X 


xAy 


^max{0,a;P + yP - 1} 


min{a;, y} 


X ■ y 


xVy 


minjl, + yP\ 


max{a;, y} 


min{l, X + y} 


x =J> y 


ii X < y then 1, else 1 — x^ + y^ 


max{l — X, y} 


min{l, 1 ~ X + y} 





Godel logic 


Product logic 


Probabilistic logic 


^x 


if a; = then 1 else 


if a; = then 1 else 


1-x 


X A y 


min{a;, y} 


X ■ y 


X ■ y 


xVy 


max{a;, y} 


X + y — X ■ y 


X + y — X ■ y 


X => y 


ii X < y then 1, else y 


ii X < y then 1, else x/y 


y - X + X ■ y 



these functions is continuous on [0, 1] or [0, 1]^ only when the underlying logic is 
neither a Godel or Product logic. The implication function on [0, 1]^ is continuous 
only for the Zadeh, CFMS or Probabilistic logic. Thus, although each of these 
6 logics is a continuous logic, only a few of them define continuous functions on 
[0, 1]^, per se. This last statement has a bearing on the complexity of cellular 
automata based on multi-valued logics. 

We begin by defining simple boolean (two state) automata based on standard 
logic. A (boolean) cellular automaton is a collection of black and white cells 
(or ones (1) and zeros (0) resp., or true (T) or false (F) statements, resp.) 
arranged on an infinite one-dimensional strip. We also call this an initial string. 
Each cell is preceded by a cell behind it and a cell in front of it; these are 
called the neighbors of the cell. We distinguish an arbitrary (but fixed) special 
element of this strip by labeling it xq (or by Xq). Returning to our band of 
initial values (or our initial configuration) we define a local rule, by a function 
g : {0, l}'^ i~> {0, 1}, on consecutive triples Xi-i,Xi,Xi+i of this band with the 
effect that ,f{x)i — g{xi-i, Xi, Xi+i) where the value f{x)i is placed directly 
below Xi (on another such parallel band) and this update of the original cell 
values is performed for all the quantities in the initial string resulting in two 
parallel bands. The procedure is then repeated to generate an infinite sequence 
of bands whose structure at infinity we wish to study. Ultimately, the idea is 
to compare and contrast the effect of the various logics used in determining the 
long term evolution of cellular automata based upon them; thus the motivation 
for calling these continuous or fuzzy cellular automata. 

We now assume standard propositional calculus. Thus, when dealing with 
(classical) boolean CA's (as in [25 ) the local rule is defined by the 8 possible 
local configurations of triples that a cell detects in its neighborhood: 

000,001,010,011,100,101,110,111 ^ro,ri,--- , rr 

where each triplet represents a local configuration of the left neighbor, the cell 
itself, and the right neighbor, and the numbers are either or 1. Letting the Vi 
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take on all possible values of O's and I's will result in the definition of 2^ = 256 
different local functions. Each one of these local functions defines a rule whose 
rule name is, by definition, the value of the sum 

7 

Rule name = ^ 2*, 

see Examples |5] |6] and [7] below. Now it is easy to see that our definition of rule 
names forces a given rule name to be an integer between and 255 and, indeed, 
every such integer defines a rule name. Observe that, by definition, the string 
Trr^Tc, . . . riro is simply the binary representation of the rule name (that is, the 
representation of that integer in base 2). 

Now, it is known that every local rule of a boolean CA is expressible as a 
disjunctive normal form (abbr. DNF), i.e., 

where dij is the j-th digit, from left to right, of the binary representation of i, and 
(resp. x^) stands for -^x (resp. x). This DNF is simply a comprehensive device 
for the writing of any one of the 255 general boolean CA's using the standard 
conjunctions and disjunctions of classical propositional logic. More precisely, 
the DNF is the form to which every statement of propositional calculus can 
be reduced, consisting of a disjunction of conjunctions each of the conjuncts of 
which is either an atomic formula or its negation. 

This representation of g is also called the canonical expression for the rule. 
For example, since 3 = l-2" + l •2^ + 0-2^ + . . ., we can write 3 = (Oil) in base 2, 
so that dsi = 0,^32 = 1,^33 1. Similarly, since 5 = l-2° + 0-2i + l-22 + . . ., we 
have 5 = (101) in base 2, and d^i = 1, (^52 = 0, c?53 = 1. In the same way we can 
show that, for example, dn = 1,^42 = 0,^43 = and dji — l,c?72 = 1,^73 = 1 
(these quantities are used in Example [5] below) . 

In addition, the rule g which sends the boolean triples 

000,001,010,011,100,101,110,111 ^ 0,1,1,0,0,1,0,1 

is a local rule. Since the effect of this rule on every such triple may be written 
as 5(1,0,1) = 1, .9(0,1,0) = 1, 5(1,0,0) = 0, 5(0,0,0) = 0, etc., this rule maps 
the initial string, say, 

...1 10 11 1... 

into (the second band in the next display, and then the third band, etc.) 

...1 10 11 1... 

1 1 1 1 ... 

1 1 1 ...&etc. 

This rule 5 is an example of a boolean CA and it is such boolean CA's (see pS] ) 
and special multi- valued/fuzzy counterparts to be defined below that interests 
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us here as an application of the general theory that is developed in [17], [18] in 
the case of an underlying natural fuzzy logic. The global dynamics (or evolution 
or asymptotics) of a rule is defined by updating the initial and all subsequent 
string values according to a repeated application of the local function to the 
neighborhood of each cell. 

Example 2 Since 110 = 2 + 2^ + 2'' + 2^ + 2^ we see that boolean rule 110 where 
110 = X]I=o ^' forces = 1 only for i — 1,2,3,5,6. Use of the disjunctive 
normal form expression above gives us 

giio(xi,a;2,X3) = Vi|r,=i Aj-^i xf\ 

= (x^ii A xi^^ A x"^^") V (sf 1 A xi^^ A Xg^^) V (x^^i A x^"" A x"^^') 

y{xf'' A xf'' A xf') V (x^^i A x*^ A 
= {x° A X2 A Xg) V {x° A X2 A a;!]) V (a;? A ^2 A xl) V (x} A a;^ A xj), 

V(a;il AxaAx") 
= (^xi A -^X2 A xs) V (^xi A a;2 A ^I'a) V (^xi A a;2 A x^) V 

V(a;i A A X3) V (xi A a;2 A ^xa) (1) 

for Xi G {0, 1}, i = 1, 2, 3, in the case of a boolean logic or x^ e [0, 1], i = 1, 2, 3, 
in the case of a multi- valued/fuzzy logic. In the sequel we will assume that ^ 
is the canonical form of EGA 110 devoid of any underlying logical structure. 

One can now ask the following natural question: Under what conditions is 
EC A 110 under a finitely-valued logic universal? The corresponding result for 
classical (two-valued logic) is believed by everyone (cf., [[3S], p. 865]), but this 
is not clear for even a trivalent logic. The point is that a general method for 
determining the universality of such an EGA under any logic is lacking. 

3 EGA 110 under continuous or fuzzy logics 

We motivate continuous or fuzzy cellular automata (or continuous EGA defined 
by them) by means of an example. In this case we note that the boolean space 
consisting of the two points and 1, that is, {0, 1}, is now replaced by an interval, 
namely, the closed interval of points between and 1, denoted by [0, 1]. 

Example 3 (Fuzzy rule 110 also known as FCA 110) Consider the continuous 
rule g defined by (/(xi, X2, X3) — X2+X3 — X2X3 — X1X2X3, for Xi G [0, 1], i = 1, 2, 3. 
This is, in fact what we call fuzzy rule 110 in 115^ . This expression for g is 
obtained by casting boolean ECA 110 (see Example^ above) under CFMS logic 
(this process is called fuzzification) . That is, in the DNF we redefine for real 
numbers a,6 e [0, 1], the quantities (a V b) to be (a + b), (a A b) to be (ab), and 
(^a) to be (1 — a), in accordance with CFMS logic. In other words, a\/ b = a + b, 
a A b = a ■ b, and -^a — 1 — a, where + and "■ " are ordinary addition and 
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multiplication of real numbers. Specifically, starting from the DNF of EGA 110, 
we find (after fuzzification) 

g(^Q{xi,X2,x:i) = (^xi A A x^) V {-^xi A X2 A ^xa) V (^xi A X2 A X3) V 
V(xi A A X3) V (xi A X2 A ^xs) 

= (1 - Xi)(l - X2)X3 + (1 - Xi)x2(l - X3) + (1 - Xi)x2X3 + 

+Xi(l - X2)X3 + XiX2(l - X3), 
= X2 + X3 - X2X3 - X1X2X3. (2) 

As an example of the long-term dynamics of this rule, we let this FCA act on 
the finite string consisting of a singleton a = 0.426 in a background of zeros (i.e., 
"a" is surrounded on both sides by zeros, see Table [1] below). Since 5(0, 0, 0) = 
the value below the central "0" of any string of the form 0,0,0 in Table [ijis zero. 
Next, (7(0, 0, .426) — .426 so this value is placed under the central zero of the 
triple 0, 0, .426. Next, ff(.426, .670, .426) = 0.689. Hence this value is placed in the 
entry under the cell value 0.670. Similarly, g(.937, .700, .510) = .518 (which goes 
under 0.7) and the process is repeated indefinitely for all possible cell values 
in every band. This generates an infinite number of bands consisting of the 
"window" of values below in Table [l] (also called the space-time diagram of the 
automaton). The basic question is, "What happens to the cell values when we 
look at arbitrarily large number of such bands?" Do the cell values converge? If 
so, to what? Is there chaos? 



.0 







0.426 0. 



...0 





.426 









426 
964 








.426 
937 
.701 






.426 
.891 





,426 
,811 
,703 



700 .510 

518 .588 



0.426 

.670 
.689 
.583 
.586 
.635 



.426 
.426 
.426 
.426 
.426 
.426 



0... 









Table 1. Evolution of a nonlinear rule with initial seed a = 0.426 (truncations 
of actual values used) 



In this sense the initial string now consists of a set of fuzzy states, that 
is a collection of arbitrary but fixed real numbers in the closed interval [0,1]. 
Inherent in this procedure is the fact that fuzzification allows one to move from 
the discrete (boolean) to the continuous (fuzzy) by extending the domain of 
definition of the rule (by modifying the underlying multi- valued logic) . We now 
turn to the process oi fuzzifying boolean rule 110 under probabilistic logic. 

Example 4 (Probabilistic rule 110 or PCA 110) For real numbers a,b ^ [0, 1], 
in the DNF for rule 110 we define the quantities (a V 6) to be {a + b — a ■ b). 
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(a A b) to be {ab), and (-la) to be (1 — a), in accordance with the probabilistic 
logic above. In other words, aV b = a + b — a ■ b, a A b = a • b, and = 1 — a, 
where + and "■ " are ordinary addition and multiplication of real numbers. As 
before, using the DNF of EC A 110 as a starting point, we find 



gi-^Q{xi,X2, xz) = {^xi A ^X2 A X3) V {^Xi A X2 A ^^3) V (^Xi A 0:2 A X3) V 
V(a;i A -1X2 A 0:3) V {xi A X2 A -'Xs) 
= (l-a;i)(l-a;2)a;3(l-^) + ^ 



where A = {1 - a;i)a;2(l - a:3)(l - B) + B, B = {1 - a;i)a;2a;3(l - C) + C, and 
C = x{l — y)z + xy{l — z) — x'^yz{l —y){l — z). Written out in full this becomes: 



(1 - a;i)(l - a;2)a;3(l - (1 - a;i)a;2(l - a;3)(l - (1 - a;i)a;2a;3(l - 

Xl(l - .T2).X3 - XlX2{l - X3) + xlx2X:i{l - .T2)(l - X3)) - 
- X2)X3 - XiX2{l - X3) + xlx2X3{l - X2){1 - X3)) - 
(1 - Xl)x2X:i{l - Xl{l - X2)X3 - XlX2{l - X3) + xlx2X3{l - 
X2)il - X3)) - Xi{l - X2)X3 - XiX2{l - X3) + xlx2X3{l - X2){1 - X3)) 

+ (1 - a;i)a;2(l - a;3)(l - (1 - a;i)a;2a;3(l - xi{l - X2)x3 - xiX2{l - X3) 
+xlx2X3{l - X2){1 - X3)) - xi(l - X2)a;3 - 0:1X2(1 - X3) 

+xlx2X3{l - X2){1 - X3)) + (1 - Xi)a;2X3(l - Xi(l - X2)X3 - 
XiX2{l - X3) + xlx2X3{l - X2)(l - X3)) + Xi(l - X2)X3 
+XiX2(l - X3) - XiX2X3(l - X2)(l - X3), 
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or, equivalently, 

X2 + X3 — X1X2X3 ~ 3 x^X2X'^ — 3 XiX2X^ + 5 

+4 X2x'^Xi — 9 x^x^xi + 4 X2X3X1 + 5 X2X^Xi + 5 x^x\x\ — 13 a:2X3xJ — 

7 iZ/ 2*^3*^1 Yx^^X^X^ ~\~ 1 1 '^3*-'' 1 ^ X^X^ Xx-^j^ H~ X-^X^ X'yX'^ 

,22,922 32 23,33,23 o242 
~\~X-yX'2 ~r ^■^2'^3 *^2*^3 *^2'^3 "^2*^3 X2^X^X\ o X^Xr^X-^ 

+6 

'^2*^3'^1 3!/-j^i2'2'2'2 ~l~ 33-2"3^2'Z']^ X '2X ^X X -^^X '2X ^ ~\~ 2 X '2X ^X 

~\~^ X'2X^X-^ X23^33^]^ H~ X-^X^X^, X-^X^X2^ ^X-^X^X^ X-^X2^X^ 

H" 3^ X 2 3 ~l~ Xq^X^X^^ 3 X^X ^X ~\~ 6 X2^X^X-^ ~\~ 3 3^2'-^3*^i '^2'^3'^1 

~|~3 X2'Z'3'2'1 X^X^X-^ 8x2'2'3X2 X2^X^X^ 43^2'-^3*^x 3^-1^X33^2 

~\~ X 2^X 3 X 2 H~ 3/ 2 X 3 X 3 X 2^X ^X 3 X '2X ^X 2 X 2^X ^X + 6 3/ 2 3 ]^ 
3 4 _ 4 4 5 _ 3 3 5 , 9 5 4 3 , 9 4 3 5 _ 

4 x-1^3^33^2 x-^x^-i^^ 2 X ^x^x ^ ~\~ x^x^x^ ~\~ 3^2*^3*^1 3 X'2'i^^'^i^ 

+3 X -^X ^X 2 33/2'^3*-^l 23/2*^3'^]^ 63^2'-^3'^j^ 63'-|^3^33^2 ~l~ 2 3- 2 -^3 1 
~|~3^2'^3'^1 33^2*-'^3'^l ~l~ 33/-|^X33^2 3/-j^3^3X2 H~ "^X-^X^X^ X-^X^-^^^ 

Now the Zoofc itp iaWe of a rule g on [0, 1]'^ is the image of the 8 boolean triples 
(0, 0, 0), (0, 0, 1), ... , (1, 1, 1) under g. Given the nature of probabihstic EGA 110 
(PGA 110) as seen above, we observe that the look up tables of probabilistic 
and FGA 110 are identical. In other words, the values of the corresponding 
fuzzification g{x, y, z) all agree on the 8 basic triples. Although more difficult to 
write down explicitly, it can be verified directly without too much effort that this 
is also the case for Product EGA 110 (see Table [2] below) , Godel EGA 110 and 
even the Lukasiewicz EGA 110 (for any p > 1). In addition, these continuous 
EGA also agree with Boolean EGA 110 when restricted to the Boolean space 
{0,1}. In fact, the reason for this equivalence of look-up tables, for any rule 
number, is that all the multi-valued logics mentioned here agree with standard 
propositional calculus when restricted to the two-point set {0, 1}. 



Boolean triples 


Boolean 110 


Fuzzy 110 


Probabilistic 110 


Product 110 


000 














001 


1 


1 


1 


1 


010 


1 


1 


1 


1 


Oil 


1 


1 


1 


1 


100 














101 


1 


1 


1 


1 


110 


1 


1 


1 


1 


111 















Table 2. The look-up tables of the elementary cellular automaton, EGA 110, 
under diverse multi-valued logics 
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4 The long term dynamics of EGA under CFMS logic 

Example 5 Since 184 = 2^ + 2'' + 2^ + 2"^ we see that the rule number 184 = 
Y^1=o ^' forces ri — 1 only for i = 3, 4, 5, 7. Use of the disjunctive normal form 
expression above gives us 

gi84:{xi,X2,X3) = Vi|r, = i Aj^^i xf', 

= (x^*^^ A x"^'' A x'^''') V {x'l^' A x^^' A x"^"') V {xf'' A x^'' A 

V(a;fi Ax^^^ Axf^'), 
= (x? A A xj) V (x} A X2 A xg) V (x} A x^ A X3) V (x} A Xj A Xg), 
= (^Xi A X2 A X3) V (xi A ^X2 A ^xa) V (xi A ^X2 A X3) V 

V(xiAx2Ax3) (3) 

= (1 - Xi)x2X3 + Xi(l - X2)(l - X3) + Xi(l - X2)X3 + X1X2X3, (4) 

= xi - X1X2 + X2X3. (5) 

Note that in "fuzzifying" the DNF ([T]), we replaced ^x by 1 — x, x V y by x + y, 
and X A 2/ in Q by their product, xy, so as to find the fuzzy form (or rule) given 
by Q or equivalently in accordance with the fuzzification process defined 
above. In this case, the local rule maps the triples of zeros and ones as follows: 

000,001,010,011, 100, 101, 110, 111 ^ 0,0,0, 1, 1, 1,0, 1. 

Thus, (boolean) rule 184 is given by ([3| above while FCA 184, given by ([5]), may 
be written as 

5i84(x, y,z) = X - xy + yz, 

for any value of {x,y,z) G [0,1]'^, (the unit cube in R'^.) From now on, "rule 
n" refers simply to the local EGA correponding to (classical) boolean cellular 
automata defined above (or in ) while "FCA n" is the CFMS-fuzzification of 
the local rule as described above. We emphasize that to each classical boolean 
CA corresponds a unique CFMS CA and vice versa. 

Example 6 Rule 18 = 2 + 2^ has the local rule 

(000,001,010,011,100,101,110,111) ^ (0,1,0,0,1,0,0,0). 
Its canonical expression being 

5i8(a;i, X2, X3) = (^xi A ^X2 A X3) V (xi A ^X2 A ^X3), 
we obtain its fuzzification as 

ffi8(a;i, X2, X3) = (1 - X2) • (xi + X3 - 2x1 • X3). (6) 
One of the EC A that motivated this paper follows. 
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Example 7 Rule 30 = 2 + 2^ + 2^ + 2'' has the local rule 

(000,001,010,011, 100, 101, 110, 111) -> (0, 1, 1, 1, 1,0,0,0). 
Its canonical expression is 

g^Q{xi,X2, xz) = (^a;iA^a;2Aa;3) V(^a;iAa;2A^a;3)V(^a;iAx2Aa;3)V(a;iA^a;2A^X3) 
so its fuzzification becomes 

9i.o{xi,X2,x^) = a;i + a;2 + a;3 - 2a;ia;2 - 0:2X3 - 2x10:3 + 2x1X2X3. (7) 

The usual fuzzification of the expression aV6 is max{l, a + 6} so as to ensure 
that the result is not larger than 1. Note, however, that taking (a + h) for the 
CA fuzzification does not lead to values greater than 1 since the sum of all the 
expressions for rule 255 is 1 (i.e., g255(x,j/,z) = 1), and so every (necessarily 
non-negative) partial sum must be bounded by 1. Since every FCA is obtained 
by adding one or more of these partial sums it follows that every FCA is bounded 
below by and above by 1. 

Indeed, we note that there are many ways of fuzzifying a rule so that this 
process of fuzzification is not unique by any means (see [2T] for other possible 
fuzzifications) . In particular, we mention one type of unusual fuzzification where 
the value a V 6 is replaced by 1 — (1 — a)(l — &). If adopted, this would lead to 
another continuous rule where considerations similar to those in |17j . |18j may 
be applied to deduce rule asymptotics. In order to distinguish this definition 
from other possible definitions in the literature (e.g., |21j ) we sometimes add the 
prefix CFMS (an acronym) in honor of its discoverers, 1^. In order to make this 
report as self-contained as possible we will define an exceptional fuzzy cellular 
automaton herewith: 

Definition S A ( CFMS CA ) or FCA, for short, with local rule g is said to be 
exceptional provided its diagonal function d, is identically zero on [0,1] (where 
d{x) = g{x,x,x) — x). That is, the set of fixed points of g is the closed interval 
[0,1]. 

A complete analysis of all 255 (CFMS) FCA's shows that there only nine 
such exceptional fuzzy CA, that is, FCA 170, 172, 184, 202, 204, 216, 226, 228, 
and 240. The definition of these nine exceptional FCA, of which six are non- 
trivial, indicate that the general methods presented in [17], [18] require major 
modifications. Of these nine, three are trivial, that is fuzzy ECA 170, 204 and 
240 have simple dynamics, for any initial string, finite or infinite. The point is 
that all but possibly these six remaining FCA have non-chaotic dynamics at 
infinity (compare [3]). These inherent difficulties in the remaining six FCA were 
treated recently in a series of separate papers [19] and [20], thus completing the 
program of the classification of the dynamics of all 255 FCA under fuzzy logic. 

To begin with, we fix the notation. We recall some definitions from |7]. The 
space-time diagram from a cell x* is the set of {x*^'' | p > Oandj € {i — 
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p, • • • , i+p}}. In the case of a single fuzzy seed, a, in a background of zeros (i.e., 
its neighbors are all zero), the space-time diagram is the infinite triangle whose 
vertex is at the singleton a = Xq. Thus, a;™„ will denote the cell at ±n steps to 
the right/left of a (the zero state) at time m. For example, x^i — 0.5429, cc^ = 
0, — 0.875, ... in the adjoining Tablejsjfor the evolution of FCA 110 starting 
from a single seed a = 0.5 in a background of zeros. The single cell Xq will always 
be denoted by a and generally we will take it that < a < 1, since a = leads 
to trivial asymptotics. 



Time 


... _3 


Local states 
-2 -1 1 2 3 • • • 


1 


... 


.5 .5 ••• 


2 


... 


.5 .75 .5 ■•• 


3 


••• .5 


.875 .6875 .5 • • • 


4 


••• .9375 


.6601 .5429 .5 ■ • ■ 



Table 3. FCA 110: Partial evolution from the point a = 0.5 in a zero back- 
ground. 



By way of further motivation for the study of the global evolution or temporal 
asymptotics of fuzzy EGA we assume that the initial string is the seed a — 1, 
in a background of zeros. For a given integer n where 1 < n < 255 it is readily 
verified that the space-time diagrams of both the fuzzy and boolean rule n are 
identical. In other words, restricting our given FCA n to the case where a = 1, 
the resulting evolution of O's and I's in its space-time diagram will necessarily 
produce an output that is identical to that which we would obtain had we applied 
(boolean) rule n to our initial string. The reason for this is that (the name of) 
a FCA is uniquely determined by its action on the eight elements 000, 001, . . . 
so one need only calculate these eight values using the fuzzy local function and 
compare them to the corresponding boolean values obtained by using the DNF. 

In addition, proceeding from a single seed in a background of zeros again, the 
space-time diagram produced by applying FCA n to this string will "converge" 
to the boolean one as a 1^. In this way we show that complexifying the 
problem by introducing a fuzzy seed (as opposed to a boolean one) can actually 
lead to a somewhat deeper understanding of the evolution of a given boolean 
rule. 

In this respect we use the following device (see [7]). Fix a FCA and choose 
an initial seed a = Xg, in a zero background such that a is very close to, but 
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less than, the number 1. We then fix a row number or, equivalently, a specified 
time T > 0. Consider the fuzzy evolution of the rule up to time T. Thus, we 
are looking at the light-cone proceeding from a to the row whose elements are 
a^^'p, x^Tj^ii ■ ■ ■ , Xq, . . . , x'^^^, xTp. Because of the results in |T7], |18j on the de- 
terminacy of the asymptotic limits of all but six FCA we have that the limits 
are continuous functions of the seed a. In addition, we recall that both boolean 
and fuzzy space-time diagrams are identical if a = 1 . We now refer to the fuzzy 
space-time diagram: If any cell at time no larger than T has an entry greater 
than 1/2 we color it black. On the other hand, if a cell at time no larger than T 
has an entry less than 1/2 we color it white. The result is that if the seed a is 
chosen sufficiently close to 1 , then the space-time diagrams of both the boolean 
and the fuzzy EGA up to time T will be identical Thus, one can expect that the 
asymptotic evolution of the fuzzy space-time diagram will have bearing on the 
distribution of O's and I's in the corresponding asymptotic boolean evolution. 

This method is illustrated in Table |4] and Table [s] below and adapted to 
the space-time diagram for FCA 110, emanating from a single seed, for sim- 
plicity. We recall that the local function for this rule is given in [TS], or that 
9iio{x,y,z) = y + z~yz~xyz, for any x,y,z G [0,1]. Thus, 5iio(.99, .09, .95) = 
.86, 5iio(.14, .20,.25) = .4, etc. 



.95 .95 
.95 .99 .95 
.95 .99 .09 .95 
.95 .99 .90 .86 .95 
.95 .99 .14 .20 .25 .95 
.95 .99 .86 .28 .40 .91 .95 
.99 .17 .65 .47 .84 .64 .95 



Table 4. Partial space-time diagram of FCA 110 from a = 0.95 up to T < 7. 



1 

11 

10 1 

1111 

10 1 
110 11 

10 10 111 



Table 5. Partial space-time diagram of boolean ECA 110 from a = 1 up to 
T <7. 
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Remark 2. Note that we truncate the values of the cells in Table S] above to two 
decimal places since we are only interested in values relative to the number 1/2. 
Once that cell value is greater than 1/2 we color it black (or give it the value 1) 
and if it is less than 1 /2 we color it white (or give it the value 0). In this way we 
see that the space-time diagrams described in Tables are identical. 



Remark 3. At this time we cannot decide the color of a fuzzy cell whose value is 
exactly one-half. Although this cell value equal to 1/2 is unlikely in the case of 
random inputs of the initial seed a, it can occur if, say, the seed value a = 1/2 
(such as in FCA 90, [7J). Indeed, it is a consequence of the definition of fuzziness 
here that given any one of the 256 fuzzy EGA, the closer the seed value a is 
to 1 the more "agreement" there is between the boolean and fuzzy space-time 
diagrams, using the device mentioned above that allows the transition from a 
fuzzy scenario to a boolean one. Consequences of this remark will be addressed 
at a later time in a different work. 



4.1 The Dynamics of FCA 110 

Interest in the boolean form of FCA 110 has been tremendous (see (5S]) since 
Matthew Cooke's announcement that it can support Turing-complete computa- 
tion. We use the general theory developed in [17], [18] to derive the dynamics 
of the fuzzification of this rule and so provide a simple proof of the results ob- 
tained in ^15j. We note that the long term dynamics of FCA 110 were obtained 
previously in |15| using special arguments pertaining to the form of the rule 
itself along with its representation as a Taylor polynomial, reminiscent of the 
techniques used earlier in [7]. 

The canonical expression for this rule is given by 

giw{x, y,z) + z-yz- xyz. 

The fixed points in U of this FCA are x — 0, repelling, x — {V5+l)/2, which lies 
outside the domain of the rule and finally, the golden number, x = (\/5 — l)/2 
which is attracting. The dynamics of this rule are given by solving the equations 

■^m+i(a) = .9110 {L:^n-lia),L-^{a),L;^_^_l{a)) 

using the initial conditions Lq = a, Lj~ = 1, (see [18] for a proof). So, all limits 
are found recursively with the first few limits being given by 

L-- ^ L--^ L--"+^ i-- " + 2 L--^ 

^2 - l + a' ^3 - 2' ^4 - ^^2' ^5 - 3' ^6 - 2a + 3' ^7 - 5, 

and 

_2a + 3 

An induction argument shows that the even and odd subscripted left-diagonal 
limits and i^„_i are given more simply by 
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for each n > 1. 



Remark 9 As an aside we note that the second formula in ^ defines the con- 
tinued fraction expansion (cf, YM^) of the golden number, cj) = (\/5 — l)/2 = 
0.618033 . . .. As a result, the fractional representation o/L^„_|_]^ appearing in ^ 
must be quotients of Fibonacci numbers (viz., i^„_i = Fn/Fn+i)- where F„ is 
the n^^ Fibonacci number, F^ = Q. It follows that the golden number is a limit 
point of the infinite sequence of left-diagonal limits in the evolution of FCA 110. 
Next, one can show by induction that 



j^— ^ F^i—\ Fji 

a Fn + Fn^i 



for each n > 1. 



Factoring out a copy of Fn from both numerator and denominator in ^ , passing 
to the limit as n ^ oo and using standard Hmiting properties of quotients of 
Fibonacci numbers we see that (j), independently of the seed, a. We can 

therefore conclude that the limit of the left-diagonal limits in FCA 110 actually 
exists and is equal to the golden number (see [15) for more details). 



Remark 10 These results about the limits of FCA 110 indicate that there is 
some "order" to the evolution of the space-time diagram of this specific CA. 
However, since each such limit is different from the preceding one this also im- 
plies that, generally speaking, the limit of a sequence of the form {x^} in its 
space-time diagram will NOT exist at infinity. There are FCA for which even 
these arbitrary limits DO exist, but they are few and shall be taken up in a future 
paper. 



5 The dynamics of PCA 110 

The most difficult investigation using any means whatsoever appears in this 
section. Here we begin with EGA 110 and change the underlying logic to a 
probabilistic one (see Example |4] for its local rule) and proceed to study the 
dynamics of its space-time diagram. First we note that the diagonal function of 
EGA 110 under probabilistic logic is given by the polynomial (see Example |4] 
above) of degree 15, 

+Qx^ + ibx^ - Ax^° + lQx'' + 2x-ix^ + I2x^ - 17 x^ . 

Its fixed points correspond to the zeros of gii^ix, x, x) — x, that is, the solutions 
of 

x-2x^ -5x^ + 12x^ -17x^ + 6x^ + I9x'^ -AOx^ + 35x^ -Ax^° - 

27x^^ + Ux^^ -21x^^ + 7x^'^ -x^^ =0. 
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By inspection we note that a; = is a root and there is an irreducible factor 
x'^ — x + 1 in the display above. It follows using the theory of Sturm sequences [[5], 
p. 198] that this diagonal function has exactly three real roots only two of which 
are within [0, 1] (including the one at x = 0) the other being at approximately 
xq — 0.4845. It is readily verified that there is a repelling fixed point at x = 
and an attracting fixed point at x — xq (the precise value of which is unknown). 

We now apply PCA 110 to the case of a single (fuzzy) seed a g (0, 1) in 



a background of zeros (see Section 4.1 above in the case of FCA 110). The 



nature of this cellular automaton based on probabilistic logic makes it such 
that gfiQ{x,y,z) introduced earlier is now an algebraic function (specifically, 
a polynomial in three variables of degree 5 in each variable separately). Still, 
it is possible to a certain extent to determine the long term behavior of the 
diagonals for this specific rule. For example, let us proceed as in [18] in 
the single seed case, where we choose and fix a seed a G (0, 1). Observe that 
the space-time diagram of PCA 110 is similar in shape to the FCA 110 case 
(no evolution to the right of the central column with apex at Xq). Since gf^Q : 
[0, 1]"^ [0, 1] is continuous on this domain, the method in [171 applies to give 
the limits recursively. Thus, noting that Lq = and = a, we find that 

L2{a) = 9 {L^ia),L^{a),L^{a)) 
= 9 (0,0,^2^(0)) 

= a + L^{a) -2aL:^{a) + 2a2(L^(a))^ - a3(i-(a))^ - o^iLl^ia))^ 
+a\L^{a))\ 

a cubic with real coefficients in the unknown (a) . For a non-trivial seed value 
(i.e., a ^ 0) our limit L = i^(a) must satisfy the cubic equation 

p{L) = [a^ - a)L^ + (2a - 0^)1^ - 2L + 1 = 0. (10) 

However, since p(0)p(l) < and a G (0,1), Bolzano's theorem guarantees for 
each a the existence of a real root xq^o) = L(a) S (0, 1). For such an a G (0, 1), 
an application of the implicit function theorem indicates that L2 is a continuous 
function of a at least locally. Varying the seed we get that L2 is continuous for 
every a € (0, 1). 

Chaotic evolution may occur if other zeros of p{L) present themselves in 
(0,1), since then it may be conceivable that seeds close to one another may 
converge to different limits under the action of PCA 110. However, we show 
presently that this is not the case, that is, for any given value of the seed a G (0,1) 
the cubic polynomial p{L) above has exactly one root in (0, 1). To this end let 
XQ{a) denote a real root of p{L). We follow the root as a varies. Since 

(a^ - a)xl{a) + (2a - a^)a;^(a) - 2a;o(a) + 1 = 0, (11) 
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the implicit function theorem gives that xq is a difFerentiable function of a and, 
differentiating both sides with respect to a, we obtain after some simphfication, 

3axl{a){a-l) + 2axo{a){2-a)-2^'^+xl{a) {{2a- l)xo{a) + 2 - 2a) = 0. 

(12) 

First note that the root cannot be stationary and, in fact, the motion of this 
root is monotone, for otherwise, we must have dxo/da — for some a which, 



along with ( 12 1 forces Xq = {2a — 2) /{2a — 1). However, since this a G (0, 1), this 
last equation is not possible for our root in (0, 1) and this establishes our result. 

Now we prove the uniqueness of the real root xo{a) S (0, 1). If possible, let 
p{L) have two real roots; then clearly, all its roots xo,xi,X2 must be real. The 



classical relationship between these roots and the coefficients of (111 gives the 
result 

2x0X1X2 = a;oa;i + a;ia;2 + a;oa;2. (13) 

An application of the arithmetic-geometric inequality {A+B+C)/3 > {ABC)^^^, 
(see [T2]) for the case of three positive numbers A — xqXi, B — X1X2, C — xqX2 



along with (13) gives the impossible result 8x0X1X2 > 27 since all Xi < 1. It 



follows that no two real roots can exist for ( 10 1 and so the real root Xo(a) found 
above is unique. The consequence of all this argumentation is that the limit 
L^(a) is well-defined, continuous and unique (it can be found numerically using 
Newton iterations) and thus no sensitive dependence (see [T]) on the initial seed 
a can occur for this diagonal. It is likely the case that the same holds true for 
the limits of all subsequent diagonals but we cannot prove this (due to the huge 
number of terms that present themselves). Nevertheless, we suspect that chaos 
cannot occur in the evolution of a space-time diagram from a single seed in 
a homogeneous background of zeros for PCA 110 although we note that this 
apparent lack of complexity seems to have no bearing on whether PCA 110 is 
universal or not (the answer to which is unknown to us). 



6 Concluding Remarks 

Now that the program of investigating the dynamics of all fuzzy cellular au- 
tomata initiated in [17], [18] is completed with the appearance of [19], [20] . 
attention is turned to studying the dynamics of these 255 EGA relative to modi- 
fication of the underlying logical structures. The so-called boolean case (induced 
by a two- valued logic) has been studied extensively beginning with the works of 
Wolfram (see e.g., [25] and the references therein). Specializing to EGA 110 we 
show that although boolean EGA 110 is universal, its fuzzy counterpart, EGA 
110, cannot be as its space-time diagram cannot be chaotic. In addition, modi- 
fying the underlying logic to a probabilistic one, as we have done in this paper, 
the techniques introduced indicate that such asymptotic data can be calculated 
at least implicitly if not explicitly, the surprising result being that even a prob- 
abilistic logic applied to EGA 110 seems to exhibit non-chaotic behavior! To 
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what extent chaos is a part of boolean CA's under ordinary logic is uncertain, 
but it appears that some CA's under continuous logics may surprisingly exhibit 
non-chaotic space-time diagrams. 
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